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I report on the effects of a periodic modulation of the con- 
trol parameter on electroconvection in the nematic liquid crys- 
tal 152. Without modulation, the primary bifurcation from 
the uniform state is a direct transition to a state of spatiotem- 
poral chaos. This state is the result of the interaction of four, 
degenerate traveling modes: right and left zig and zag rolls. 
Periodic modulations of the driving voltage at approximately 
twice the traveling frequency are used. For a large enough 
modulation amplitude, standing waves that consist of only 
zig or zag rolls are stabilized. The standing waves exhibit 
regular behavior in space and time. Therefore, modulation 
of the control parameter represents a method of eliminating 
spatiotemporal chaos. As the modulation frequency is varied 
away from twice the traveling frequency, standing waves that 
are a superposition of zig and zag rolls, i.e. standing rect- 
angles, are observed. These results are compared with exist- 
ing predictions based on coupled complex Ginzburg-Landau 
equations. 

47.54.+r,05.45.Jn 



I. INTRODUCTION 

When a spatially extended system is driven far from 
equilibrium, a series of transitions occurs as a function 
of the external driving force, or control parameter [0. 
The initial transition is typically from a spatially uniform 
state to a state with periodic spatial variations, called a 
pattern. As the control parameter is increased, a se- 
quence of instabilities occurs that produces increasingly 
complex spatiotemporal patterns. Ultimately, the system 
becomes fully turbulent. States of spatiotemporal chaos 
form an interesting class of patterns |||,||] . Loosely speak- 
ing, spatiotemporal chaos refers to deterministic patterns 
that possess a random variation in space and time. They 
generally exhibit an underlying periodicity but are char- 
acterized by a finite correlation length and correlation 
time. Since their discovery in fiuid dynamical systems 
[^, states of spatiotemporal chaos have been observed 
in a wide range of systems, including Rayleigh-Benard 
convection, Faraday instabilities, Taylor- Couette fiow, 
electroconvection in nematic liquid crystals, lasers, and 
chemical reactions Despite the ubiquitous nature 

of the phenomenon, a full understanding of spatiotem- 
poral chaos remains one of the outstanding problems in 



nonlinear dynamics. 

One of the challenges facing the study of spatiotem- 
poral chaos is the difficulty associated with characteriz- 
ing the dynamics [|l|. Because the systems are spatially 
extended, they are inherently high dimensional. This 
precludes the use of many of the successful techniques 
developed to study low dimensional chaos in dynamical 
systems [Q. Also, most examples of spatiotemporal chaos 
occur at sufficiently large values of the control parame- 
ter that weakly nonlinear techniques, such as amplitude 
equations, are only applicable as phenomological models, 
if at all. The work reported here uses a state of spa- 
tiotemporal chaos that occurs in electroconvection in 152 
[||. In contrast to other systems, this state does occur in 
the weakly nonlinear regime where quantitative compar- 
ison between theory and experiment is possible. A fun- 
damental theory of electroconvection, the WEM model 
[^,0 , exists and provides the necessary starting point for 
a quantitative derivation of the relevant amplitude equa- 
tions Therefore, this system is an ideal candidate for 
studying spatiotemporal chaos. However, currently only 
two of the four necessary amplitude equations have been 
derived and only qualitative predictions exist In this 
work, I use temporal modulation as a probe of the sys- 
tem's dynamics. This provides both a test of the existing 
amplitude equation description and a means of guiding 
future experiments and calculations. 

For electroconvection |p|,p!Q[, a nematic liquid crystal 
is placed between two properly treated glass plates so 
that the director is everywhere parallel to the plates and 
along a chosen axis. A nematic liquid crystal is a ffuid in 
which the molecules have orientational order, and the di- 
rector refers to the axis parallel to the average alignment 
of the molecules [Q. The liquid crystal is doped with 
ionic impurities, and an ac voltage is applied across the 
sample using transparent electrodes on the glass plates. 
Above a critical value of the applied voltage, there ex- 
ists a transition from uniform conduction to convection 
rolls, with a corresponding periodic variation of the di- 
rector and concentration of ionic impurities. Because the 
system is anisotropic, the patterns can be characterized 
by the angle 9 between the director and the wavevector 
of the pattern. Electroconvection in the nematic liquid 
crystal 152 exhibits a forward, Hopf bifurcation to oblique 
rolls 1^ . A Hopf bifurcation is a transition to a traveling 
wave pattern, and oblique rolls correspond to patterns 
where 0° < ^ < 90° loll. At onset, because the director 



1 



only defines an axis and not a positive direction, states 
with wavevectors of equal magnitude but with angles 6 
and TT — 6 are degenerate and referred to as zig and zag 
rolls, respectively. The state close to onset that is stud- 
ied here consists of four modes: right- and left-traveling 
zig and zag rolls ||]. It is the interaction of these four 
modes that results in spatiotemporal chaos. 

Resonant modulation of the control parameter in a sys- 
tem with a Hopf bifurcation is known to stabilize stand- 
ing waves for large enough modulation strength [|l2| -p!7[. 
However, it has not been applied previously to a state of 
spatiotemporal chaos. There are three main questions 
addressed in this work. First, I have made a survey 
of the range of existence of standing wave patterns and 
their qualitative features for a reasonable set of values of 
the applied voltage, applied frequency, the modulation 
strength, and the modulation frequency. Second, there 
are three possible standing wave solutions: standing rolls 
(only zig or zag rolls are present), standing rectangles (su- 
perimposed zig and zag rolls with equal amplitude), and 
standing cross rolls (zig and zag rolls with different am- 
plitudes) I will show that the standing wave states 
are generally standing rolls, but that standing rectangles 
can be observed. Finally, I will demonstrate that the 
standing roll states are temporally regular and that they 
eventually become spatially uniform. Therefore, this is 
an example of the elimination of spatiotemporal chaos 
and provides an interesting contrast to systems where 
temporal modulation produces irregular behavior in an 
otherwise regular system |ll5| , p!8[ . 

The rest of the paper is organized as follows. Section II 
provides the details of the experimental techniques. Sec- 
tion III presents the experimental results. In this section, 
I will report separately on the results of applying mod- 
ulations of the control parameter below and above the 
critical voltage for the onset of convection in the absence 
of modulations. Section IV will discuss the relationship 
between these results and existing predictions of relevant 
amplitude equations. 



II. EXPERIMENTAL DETAILS 

The experiments were carried out using two electro- 
convection cells containing the liquid crystal 152 [ p!9| . 
The first cell was a custom made cell with a thickness 
of 25 fim. It was formed from two glass slides that were 
coated with a layer of indium-tin oxide (ITO), a trans- 
parent conductor. The conductive coating was etched to 
form a 0.5 cm x 0.5 cm square electrode in the center of 
a 2.5 cm x 2.5 cm cell. The director was aligned using a 
rubbed polyimide. The second cell was a commercial cell 
obtained from EHC, LTD. ||2^]. This cell had a thickness 
of 23 jam. The electrode was a 1.0 cm x 1.0 cm area of 
ITO in a 2.5 cm x 2.5 cm cell. The alignment was also 
due to a rubbed polyimide coating. For both cells, there 
was some forcing of a pattern due to fringing fields at 



the edge of the electrodes. However, a region existed in 
the middle of each cell where convection started sponta- 
neously. All measurements were made in this region of 
the cell to minimize boundary effects. 

The main differences between the two samples were the 
values of the critical voltage and the Hopf frequency due 
to differences in iodine doping. The custom sample was 
filled with 152 that had been doped with 6.3% by weight 
molecular I2 two months prior to filling the cell. The 
commercial cell was filled with 152 that had been doped 
with 6% by weight molecular I2 seventeen months prior 
to filling. The custom cell was aged for five months after 
filling before experiments were started, and the commer- 
cial cell was aged for one month after filling. For the 
custom cell, the critical voltage was approximately 21 
V and the Hopf frequency was approximately 0.125 Hz 
at an applied frequency of 25 Hz. For the commercial 
cell, the critical voltage was approximately 15 V and the 
Hopf frequency was approximately 0.34 Hz at an applied 
frequency of 25 Hz. 

It is known for electroconvection in 152 that the level 
of iodine in the sample will drift in time. This re- 
sults in drifts in the critical voltage and the Hopf fre- 
quency. Because of this drift, the following protocols 
were used. The critical voltage and Hopf frequency 
were measured before and after every set of measure- 
ments. The modulated voltage had the following form: 
V{t) = [Vo + VmCos{umt)]cos{rtt) . The two main di- 
mensionless control parameters are: e = /V^ — 1 and 
b = Vm/Vc. Here Vc is the critical voltage for the onset 
of convection in the absence of modulation. The drift in 
Vc was linear in time and corresponded to a drift in e 
of 0.001/Hr. This drift is accounted for in all reported 
values of b and e. The modulation frequency will be dis- 
cussed in terms of the shift from resonance (/* — /m/2), 
where fm is the modulation frequency and /* is the nat- 
ural frequency of the pattern. For e < 0, /* is the Hopf 
frequency and for e > 0, it is the frequency of the 
pattern in the absence of modulation. Despite the dif- 
ferences between the two cells and the drift in time of 
Vc, the behavior as a function of 6, e, and (/* — frnf^) 
is completely reproducible. The sample temperature was 
held constant at either 40 =b 0.002 °C or 42 ± 0.002 °C. 
The later temperature was used after the total drift in 
Vc at 40° C had exceeded approximately one volt. 

Images were taken using a standard shadowgraph 
method and are presented here with the undistorted 
director aligned in the horizontal direction. Because 
of the well-known nonlinear effects of the shadowgraph 
method Q , the images have been Fourier filtered so that 
only the fundamental modes are present. This is essen- 
tial for highlighting the standing wave character of the 
modulated pattern. 

The frequency of the pattern was determined by tak- 
ing the Fourier transform of a time series of 32 images. 
The images typically covered a spatial area containing 
18 rolls, though smaller regions of only 7 rolls were also 
used. The time between images was chosen so that the 
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time series covered 4 to 6 periods of the fundamental 
frequency. 

For the measurements of the dynamics of the local am- 
plitude of each mode, a time series of 32 images was used. 
The time between images was chosen so that the series 
consisted of four cycles of the pattern. Each image cov- 
ered a spatial area of approximately 7 wavelengths. The 
modulus squared of the space-time Fourier transform of 
the series, 5'(k, cj), was used to compute the amplitude 
of each mode. The power in each mode, right- and left- 
traveling zig and zag rolls, was determined by summing 
5'(k, uj) over a 5 by 5 pixel grid in wavenumber space and 
5 pixel window in frequency space. The grid and window 
were centered on the peak in 5'(k, uj) that corresponded 
to the mode of interest. The amplitude of each mode is 
the square root of the power. 

For measurements of the onset of standing waves, e 
was fixed and the value of b was either stepped up or 
down in increments of 0.005. At each step, the system 
was equilibrated for 10 minutes before a time series of 
images was taken. The time series of images were used 
to determine if the pattern was frequency locked to the 
modulation and whether or not a standing wave had been 
established. 



III. EXPERIMENTAL RESULTS 




FIG. 1. Four images of the pattern in a 0.6 mm x 0.6 mm 
region of the cell. The bar in (c) represent 0.2 mm. (a) An 
image of the cell at e = 0.03 and no modulation. The image 
has been Fourier filtered so that only the fundamental mode 
remains. Image (b) through (d) are three images from a time 
series taken at e = 0.036 and b — 0.04. The images are 0.9 s 
apart and have also been Fourier filtered. The modulation 
frequency was 0.694 Hz, which corresponds to twice the Hopf 
frequency. These images illustrate the standing wave nature 
of the pattern. 

Figure 1 is a comparison of the typical pattern in the 
region where spatiotemporal chaos exists and the stand- 
ing wave pattern that is stabilized by the modulation. 
Figure la is a single snapshot of the state of spatiotem- 



poral chaos at a value of e = 0.03 and no modulation. 
Figures Ib-d are three images taken 0.9 s apart of the 
standing wave state at an e = 0.036, b = 0.04, and 
fm = '^fh' The time between the images was chosen 
to highlight the relative change in phase that is charac- 
teristic of a standing wave as one crosses the minimum in 
intensity. For example, the intensities in the lower right 
corner in Fig. lb are opposite those in Fig. Id. 
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FIG. 2. The symbols give the location of the transition 
to standing waves as a function of modulation strength b and 
reduced control parameter e. The solid symbols and x's are 
for increasing the value of b. The open symbols and the +'s 
are for decreasing the value of b. The difference between the 
symbols is described in the text. For e < 0, there is no pattern 
present below the symbols, and a spatially uniform pattern 
exists above the symbols. For e > 0, the state of spatiotem- 
poral chaos exists below the symbols. Above the symbols, 
the system is phase locked to the modulation, and a state 
of standing rolls exists. For negative e, the solid line is the 
expected onset based on coupled complex Ginzburg-Landau 
equations, 6 = e. For positive e, the solid line is the curve 
b = 0.57e — 0.006 and provides a guide to the eye for the 
transition values. 

Figure 2 summarizes the range of existence of the 
standing wave patterns when a modulation of twice the 
Hopf frequency is used. The solid symbols and x's are 
the location of the transition to standing waves as mea- 
sured by increasing e. The open symbols and the +'s are 
the location of the transition to standing waves as mea- 
sured by decreasing e. In the region above and between 
the two solid lines, the pattern is composed of uniform 
standing rolls. As a check on the stability of the standing 
rolls, two runs were made at fixed b. For these runs, e 
was increased in steps of 0.005, with a waiting time of 
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5 minutes per step. One run was at b = 0.02, and the 
other was at b = 0.05. For the entire range of e within the 
boundaries, only uniform standing roUs were observed at 
these two values of b. 

For negative values of e, the behavior of the system 
is straightforward. The system makes a transition di- 
rectly from a uniform state to a state of frequency locked 
standing waves. The pattern consists of either standing 
zig rolls or standing zag rolls, and never the superposi- 
tion. The pattern is also frequency locked to half the 
modulation frequency. As shown in Fig. 2, within the 
resolution used here, there is essentially no hysteresis in 
the transition. 



less well-defined. For the purposes of Fig. 2, the onset 
was taken to be the point where the patches of standing 
rolls had a size on the order of 18 wavelengths. 
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FIG. 3. Transition to standing waves for e = —0.054. 
Below the curve, the system is uniform. Above the curve, the 
state of the system is standing rolls. The symbols are the 
experimental values. The solid line is a fit to the expected 
curve: = + (2nTdf - fm/^f • In this case, /* = 0.162, 
which is the Hopf frequency of the system. 

Figure 3 shows the behavior of the system when the 
modulation frequency is varied away from twice the Hopf 
frequency at negative e. The solid squares represent the 
onset to standing waves, and the solid line is a fit to a 
parabola. In this case, the parabola is centered on the 
Hopf frequency. 

For positive values of e, the situation is more compli- 
cated because the ground state is the state of spatiotem- 
poral chaos. In this case, I have measured the transition 
in two different ways. First, I have considered the onset 
to standing rolls. These are the solid and open squares 
in Fig. 2. When stepping down in e, the transition from 
standing rolls to the disordered state was easily identi- 
fied. However, because of the spatial disorder, when e is 
increased, the onset to spatially uniform standing rolls is 
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FIG. 4. Measured frequency of the pattern as a function 
of e. Illustrates the decrease in the pattern frequency with 
increasing e. 

A better measure of the transition for positive e is to 
use the point at which the pattern becomes frequency 
locked to the modulation frequency. This boundary is 
given by the x's and the +'s in Fig. 2. Figure 4 shows 
the plot of the pattern frequency as a function of e in the 
absence of modulation. From this, one can see that it is 
easy to distinguished the locked and unlocked patterns 
for e > 0.01, as the frequency of the pattern differs sig- 
nificantly from the Hopf frequency. With this definition 
of the transition, there is no measurable hysteresis. 

Figure 5 shows the effect of varying the modulation 
frequency for positive e. As with negative e, the solid line 
is a fit to a parabola. In this case, the parabola is centered 
on the frequency of the unmodulated pattern, not the 
Hopf frequency. Also, when (/* — /m/2) > 0.04, the 
standing wave pattern at onset is standing squares. This 
is illustrated in Fig. 6. For these images, the modulation 
frequency was 0.472 Hz, and the unmodulated pattern 
had a frequency of 0.305 Hz. In this range of (/* — 
/rn/2), increasing b leads to a secondary transition from 
the standing rectangles to standing rolls. 
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FIG. 5. Transition to standing waves for e = 0.03. Below 
the curve, the system exhibits spatiotemporal chaos. Above 
the curve, the system is phase locked to the modulation fre- 
quency. The symbols are the experimental values. The solid 
line is a fit to: 6^ = a + — /m/2)^. In this case, 
/* = 0.305 Hz^ which is the frequency of the pattern at 
e = 0.03. In contrast, the Hopf frequency is 0.34 Hz. For 
I/* — /m/2| < 0.04, the transition is directly to standing rolls. 
For (/* — /m/2) > 0.04, the transition is initially to stand- 
ing rectangles, and there is a secondary transition to standing 
rolls. 




FIG. 6. Three images from a time series taken at 
e = 0.036, h = 0.04, and fm = 0.472. This corresponds to 
(/* — fm/'^) = 0.069 in Fig. 5. The images cover a region 
that is 0.6 mm x 0.6 mm and are 0.9 s apart. The bar in 
(b) represents 0.2 mm. The images have been Fourier filtered 
so that only the fundamental mode is present. These images 
illustrate the standing wave nature of the pattern. 

I made a qualitative survey of the behavior as a func- 
tion of the driving frequency. By increasing the driv- 
ing frequency, one decreases the angle between the 
wavevector of the pattern and the undistorted director 
orientation. For our system, the same qualitative be- 
havior was observed for applied frequencies up to 80 Hz. 
For 80 Hz, 6 — 10°. For modulation at twice the Hopf 
frequency, a standing roll pattern is observed. A more 
detailed study of the effects of varying d will be the sub- 
ject of future work. 

The local temporal behavior of the standing roll state 
is extremely regular. This is shown in Fig. 7. Each plot 



is a time series of the local amplitude. The amplitude 
is measured every 2 minutes, with the initial point of 
the time series taken 10 minutes after the modulation is 
applied. Figure 7a is a plot of the local amplitude as a 
function of time for e = 0.01 and no modulation. The am- 
plitudes of the right- traveling zig rolls, left- traveling zig 
rolls, and left-traveling zag rolls have been shifted from 
their true values by 0.015, 0.01, and -0.005, respectively. 
These shifts clarify the anti-correlations present between 
the various modes and the irregular variation in time. 
This behavior has been reported previously Figure 
7b is the local amplitude for e = 0.01, h = 0.02, and 
fm = '^fh- This figure shows both the regular temporal 
behavior and the establishment of standing rolls (the zag 
amplitude has gone to zero). 
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FIG. 7. (a) Plot of the amplitudes of the right-traveling 
zig rolls (x), left-traveling zig rolls (□), right-traveling zag 
rolls (o), and left-traveling zag rolls (A). The system was at 
e = 0.01 and b = 0%. The amplitudes have been shifted 
as described in the text. (b) Plot of the amplitudes of 
the right-traveling zig rolls (x), left-traveling zig rolls (□), 
right-traveling zag rolls (o), and left-traveling zag rolls (A). 
The system was at e = 0.01 and b = 2.0%. 

The development of the local, temporal order, gener- 
ally occurred in under a few minutes. In contrast, the 
spatial ordering involves extremely long time scales. It 
can take up to two hours for the standing roll domains to 
reach sizes on order of the system size. However, upon 
removal of the modulation, the disorder develops in a 
few minutes. This difference in time scales is reflected in 
Fig. 2. One sees that there is essentially no difference 
between the transition to a frequency locked state and 
the transition to standing waves measured by decreasing 
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e. However, standing rolls of a particular size occurred 
at values of b slightly above the transition defined by fre- 
quency locking. This is easily understood in terms of the 
10 minute waiting time used when stepping b. Clearly, 
the details of the spatial ordering and the multiple time 
scales involved is an interesting problem. However, it is 
outside the scope of this paper and will be the subject of 
future work. 



IV. DISCUSSION 

The transitions for negative e can be directly com- 
pared with predictions of the relevant coupled amplitude 
equations |17|. I find excellent agreement between the 
measured onset of standing waves and the predictions of 
Ref. 0. The onset to standing waves should occur when 
b = . Here the real part of fi is just e and the imaginary 
part of /i is 27rrrf(/* - /,n/2), where = -fi£ /{n'^Kii) 
is the director relaxation time. Here 71 is a rotational 
viscosity and Ku is the splay elastic constant of the di- 
rector. For a modulation frequency equal to twice the 
Hopf frequency, the onset is given by the line b = e. This 
is the solid line shown in Fig. 1. For the more general 
case, one has 6^ = + {2iTTdY{fh — fm/'^Y- The solid 
line in Fig. 3 is a fit to this equation. The resulting 
values are: e = -0.058 ± 0.005, = 0.229 ± 0.004 s, 
and fh = 0.162 ± 0.01 Hz. For comparison, the mea- 
sured values of these parameters are: e = —0.054 ±0.002, 
Td = 0.18 ± 0.05 s, and fh = 0.164 ± 0.005 Hz. 

The nature of the standing wave pattern depends on 
the nonlinear coefficients in the amplitude equations. 
The fact that I observe standing rolls at onset has im- 
portant consequences. First, the coupling coefficient be- 
tween zig and zag rolls traveling in the same direction 
has been calculated [|[. Based on this calculation, it is 
likely that standing rectangles are the stable state for the 
parameter range in my experiments. However, the con- 
dition for the stability of standing rolls does involve all 
of the nonlinear coefficients jl^], and standing rolls are 
not ruled out by the calculations of Ref. H]. Therefore, 
these experiments highlight the need for a determina- 
tion of all of the nonlinear coefficients before quantitative 
comparisons between amplitude equations and the exper- 
iments are possible. On the other hand, in the absence 
of theoretical calculations, the temporal modulation ex- 
periments provide a means to determine the coefficients 
experimentally. 

The results for positive e are in qualitative agreement 
with the predictions of Ref. jl^]. The critical value of 
b for the transition to standing waves is linear in e for 
fixed fm- However, I find b = 0.577e — 0.006 (the solid 
line in Fig. 2), and not b = e. This is not surprising 
given that the ground state of the experimental system 
is a state of spatiotemporal chaos. This pattern can only 
be described by amplitude equations that include spatial 
derivatives, and these are not included in Ref. [|l^ 



The other qualitative agreement with Ref. [|T^ is the 
behavior as a function of modulation frequency. The crit- 
ical values of 6^ are quadratic in (/* — /^/2) for fixed e. 
However, it is clear from Fig. 5 that /* is the frequency of 
the unmodulated pattern for the fixed value of e, and not 
the Hopf frequency. This is due to the shift in frequency 
with e that is illustrated in Fig. 4. 

An additional feature of the behavior at positive e that 
requires a theoretical explanation is the regular dynam- 
ics of the standing wave state. Though an incomplete 
description, the existing amplitude equation calculations 
suggest that the unmodulated state is Benjamin- Feir un- 
stable for all wavenumbers p. This provides a possi- 
ble explanation for the spatiotemporal chaos at onset. 
From the fact that the modulated state exhibits regu- 
lar dynamics, one can infer that the standing rolls are 
Benjamin- Feir stable. This situation is the opposite of 
that previously observed in electroconvection in a dif- 
ferent nematic liquid crystal [|l5[. In that system, the 
unmodulated state was stable. For high enough modu- 
lation, the modulated state was Benjamin- Feir unstable 
and resulted in irregular dynamics [|l5|. The behavior in 
that case agreed well with calculations based on ampli- 
tude equations that included the spatial derivatives [ p2| . 
A similar calculation is needed for the system reported on 
here. In particular, it will be important to determine if 
temporal modulation is a general method for eliminating 
spatiotemporal chaos, or if it is specific to this system. 
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